Travelling waves in a drifting flux lattice 
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Starting from the time-dependent Ginzburg-Landau (TDGL) equations for a type II superconduc- 
tor, we derive the equations of motion for the displacement field of a moving vortex lattice ignoring 
pinning and inertia. We show that it is linearly stable and, surprisingly, that it supports wavelike 
long-wavelength excitations arising not from inertia or elasticity but from the strain-dependent mo- 
bility of the moving lattice. It should be possible to image these waves, whose speeds are a few 
pm/s, using fast scanning tunnelling microscopy. 



PACS numbers: 74.60. Ge, 74.20. De 

It was shown in , on general symmetry grounds, that 
an ordered array of particles moving through a dissipative 
medium (e.g., a steadily sedimenting colloidal crystal or 
a flux-point lattice drifting through a type II supercon- 
ductor) is governed by dynamical equations qualitatively 
different from those for a lattice at thermal equilibrium. 
Even the long-wavelength dynamical stability of such a 
drifting lattice was shown to rest not on its elasticity but 
on the signs of certain phenomenological parameters [see 
eqns. (0) and (|) below] governing the dependence of 
the local mobility on the lattice strain. A microscopic 
calculation |Q showed (see JIJ]) that for a sedimenting 
colloidal crystal the signs were such as to lead to an in- 
stability. We know of no analogous calculation for driven 
flux lattices. In this Letter we ask: are drifting flux lat- 
tices stable? We answer this question of fundamental 
importance starting from a time-dependent Ginzburg- 
Landau (TDGL) treatment without quenched disorder. 
We find that the moving lattice state is stable, with small- 
amplitude, long-wavelength disturbances propagating as 
underdamped waves whose speed, we emphasise, is deter- 
mined by the strain -dependent mobility and the imposed 
current, and not by inertia and flux-lattice elasticity. We 
calculate the wavespeed (a few /im/s) in terms of inde- 
pendently measurable parameters arising in the TDGL 
equations. 

We begin by summarising the derivation of the coarse- 
grained dynamical equations for a drifting lattice p| and 
defining the quantities we are going to calculate. Con- 
sider a slab of type II superconductor of thickness much 
larger than the magnetic penetration depth Xh, lying in 
the xy plane, threaded by a flux lattice (spacing Xh ) 
with magnetic field along the z direction. An applied spa- 
tially uniform transport current density Jt = J(X, gives 
a Lorentz force —Jt4>y/c per unit length on a vortex car- 
rying flux 4>z, c being the speed of light. The perfect 
flux-point lattice will then acquire a constant, spatially 
uniform drift speed Vt, = MJ t cj)/c. Here M, the macro- 
scopic mobility of the lattice, is determined by dissipative 
processes in the normal core as well as by the relaxation 
of the electromagnetic and order-parameter fields in the 



region between the vortices. Any perturbation of the per- 
fect moving lattice will result in inhomogcneities in the lo- 
cal electromagnetic and order parameter fields, and thus 
to a spatially varying flux-point velocity. The mobility is 
thus a tensor which depends on the local state of distor- 
tion of the flux lattice. For a lattice drifting along — y, 
ignoring Hall effects, pinning, inertia ||, and the effects 
of lattice periodicity, the displacement field u = (u x , u y ) 
as a function of position r and time t, defined with re- 
spect to a perfectly ordered crystal, in a frame co-moving 
on average with the flux lattice, must then obey 

d t u x = vid y u x + v 2 d x u v + D T V 2 u x + D L d 2 x u x + 

D L d x d v Uy + 0(WuWu); (1) 

d t u y = v 3 d x u x + VidyUy + D T V 2 u y + D LdyUy + 

D L d x d y u x + 0(Vu\7u), (2) 

where the terms containing the phenomenological coef- 
ficients 0j Vi tx vl arise from the "hydrodynamic" in- 
teraction of the moving vortices, Dr, = M(X + 2/z), and 
Dt = M/i, A and /i being the Lame coefficients of the flux 
lattice [^) . These equations are constructed using general 
symmetry arguments and hold for any steadily drifting 
lattice at large length scales (^> Xh, for a flux lattice). 
In this Letter we calculate the coefficients Vi for the spe- 
cific case of a drifting flux lattice, from a time-dependent 
Ginzburg-Landau (TDGL) description to which we turn 
next. The importance of the {vi} for the long-wavelength 
behaviour of the drifting flux lattice is clear: V2V3 > 
yields a wavelike dispersion whereas V2V3 < a linear 
instability. 

Scaling lengths by Xh, energies by the condensation 
energy E c in a volume , the order parameter by 
its bulk mean-field value in the superconducting phase, 
times by h/E c , the magnetic field H by \/2H c where H c 
is the thermodynamic critical field, the total electrochem- 
ical potential by E c /e* where e* = 2e is the charge of the 
Cooper pair, and defining the Ginzburg-Landau param- 
eter k = Xh/S, where £ is the bare coherence length, we 
obtain the dimensionless TDGL equations (6[|| for the 
dynamics of the superconducting order parameter ^(r, t): 
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where the phenomenological kinetic coefficient F is in 
general complex, with real and imaginary parts Ti and 
T 2 respectively. 

The equation of motion for the vector potential is given 
by Ampere's law, 



V x V x A = J,, 



(4) 



where the normal and super currents are, respectively, 



ZKl 



-df,A 



(5) 



<j being the normal-state conductivity tensor. 

We work in the large k limit, where a phase-only ap- 
proximation of the TDGL equations (||) applies, and, 
for simplicity, we set the normal state Hall conductivity 
&xy — 0. We begin by writing ip in terms of an amplitude 
/ and a phase \: 

^(r, t) = /(r, i)exp[ix(r, t)}. (6) 

In terms of the gauge-invariant vector and scalar poten- 
tials, Q = A Vx/k and P = $ + <9 t x, the magnetic and 
electric fields are then, respectively, 



h = V x Q , 
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(7) 



For large k, deep in the superconducting phase, the 
amplitude relaxes rapidly to a value determined by the 
phase. We can thus solve for / in terms of \ from (|^) 
yielding, to leading order in the effective phase-only 
TDGL equation 



d tX + $ = P 
with 7i = Re [r _1 ] and 

V x V x Q = a. 
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Lastly, charge conservation — V.(J r , 
(H) and (@) leads to 

V.(cr.E) - V.Q = 

which, with (||) and (§) implies 
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The {vi} in p), which encode the change in the 
mobility of a region of the flux lattice when it is com- 
pressed or tilted, arise primarily from electromagnetic 



field disturbances, screened on the scale A# ||. Ideally, 
therefore, we should calculate the mobility of distorted 
regions on a scale Xh- However, our main concern is the 
signs of the {i>i}, i.e., in the direction of drift of a tilted 
region and in whether a denser region drifts faster or 
slower than a rarer region. To this end, we take the sim- 
plest compressions/rarefactions and tilts, namely, those 
taking place at the level of a pair of particles. This should 
give a qualitatively correct assessment of the stability and 
a reasonable estimate of the wavespeed. Indeed, our cal- 
culation shows that the ViS decrease by a factor of 10 as 
the flux-lattice spacing varies from .25 Xh to Xh, justify- 
ing post facto this nearest neighbour approximation. We 
work, therefore, with a pair of flux points moving rigidly 
with a velocity v/,, as a function of their fixed separa- 
tion vector a. For such rigid motion, time-derivatives 
can be replaced by — v^.V. Expanding (||), (^|), and jll| ) 
in powers of vl, we obtain at 0(1) the equilibrium, time- 
independent Ginzburg-Landau equations, and at O(v^) 
a set of linear inhomogenous differential equations. 

Exploiting jlO],U the invariance of the time- 
independent Ginzburg-Landau equations under an arbi- 
trary virtual displacement d, the requirement of compat- 
ibility between vx and the imposed transport current Jt 
leads, for large k and within the phase only approxima- 
tion, to the "solvability condition" for the inhomogeneous 
0(vl) equations: 

i J dS.(3 lsXd - JdXi) = 7i / (XdP)dr; (12) 

the integral on the left-hand side is over the boundary 
of the sample, = d.VJo, Xd = d.Vxcb Jo an d Xo be- 
ing respectively the supercurrent and phase field at equi- 
librium, and the subscripts 0,1 denoting the 0(1) and 
O(vl) parts respectively of the term in question. Eq. 
( |l2| ) will yield the relation between Jt = Jtx and v£, i.e. 
the vortex equation of motion. 
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FIG. 1. Coordinate system (r, 9) for two- vortex case 
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Consider a pair of identical unit vortices, in a geometry 
denned in Fig.[l] (in cylindrical polar coordinates (r, 6, z)). 
ip and <j) are the angles made by J f with a and the vir- 
tual displacement d respectively, and 9q that between 
and the negative y-axis. We assume the flux lines to be 
parallel to the z-axis and ignore the effects of line wan- 
dering. We also define cylindrical coordinates [r\,6\,z) 
and (7*2, 62, z) with their origins at the two vortices. The 
surface integral on the left-hand side of (|l2|) can be ex- 
pressed in terms of the applied transport current. At the 
boundaries the fields are effectively those of a single vor- 
tex at the origin, with twice the winding number. There- 
fore 3is(r = oo,0) = J t , Jd.e r = 2dsm(8 — </>)/(kx 2 ), 
Xd = 2d.V6> = — 2dsin(9 — <f))/r and xi = nJtr cos 9. 
Substituting these expressions into the left-hand side of 
(|l2|), and performing the angular integration we find 
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Evaluation of the right-hand side of ( |l2| ) requires solving 
for -P(r,i) from (O) which, at 0(«l), is simply 
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V 2 P-7iP = 



(14) 



Near the centre of each vortex P w — v^.Vx and \ is 
equal to the angular variable 8\ or 82 around that vor- 
tex. Therefore P w cos^i — 9a) /ri as r\ — > and 
P « u_l cos(# 2 — 9o)/r2 as r 2 — ► 0. The solution to 
eqn. (|l4|) for the vortex pair with these boundary con- 
ditions is 

P(r) = i[K 1 (ar 1 )cos(9 1 - 9 Q ) + Xi(w 2 ) cos(9 2 - 9 Q )} 

(15) 



where v — v^a and a = K^J^\ju X x- Also, 
Xd = d.Vx = d[sm((j) - 9 1 )/r 1 + sin(^ - 



2 )/r 3 ]. (16) 



Using ( |l5| ) and (16) on the right hand side of (|l2|), and 
noting that d is arbitrary, we obtain the vortex-pair equa- 
tion of motion in the form 

2-7T 

2 — (J t x z) = Av L + B(v L x z) + C(v L .a)a 

K 

-Da.(v L x z)a (17) 



in the N term: for ip 7^ or tt/2 the centre-of-mass ve- 
locity is not parallel to the driving force. 

Now consider a steadily drifting undistorted flux-point 
lattice, and focus on a nearest neighbour pair of flux 
points with initial separation vector ao. Perturb it 
slightly: a — * a + 8a, thus causing a velocity perturba- 
tion Sv. Then we can extract the {vi} by differentiating 
our two- vortex result ( [L8| ) as follows: if ao||y and <5a||x, 
then ■jpr't- = v±; if a ||x and <5a||y, then -^fj^ = v 2 \ if 

<5r„ 
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a ||x and 5a||x, then 



Sa/aQ 

V3; if ao||y and <5a||y, then 



Sa / ag = W4; We have assumed, as justified early in the 
paper, that changes in vortex velocity due to a local dis- 
tortion are local. We find: (a) V2 > 0, so that a lattice 
moving in the — y direction will veer to the right (left) 
if the horizontal crystal planes are tilted up to the right 
(left); and (b) V3 > 0, so that a local ^-compression of the 
lattice increases the velocity in the direction of the force. 
From (0) and (§), this means the moving flux-lattice is 
stable. Also, v-y > and V4 < for the coefficients con- 
trolling the wavespeeds along the drift |llj . The resulting 
mode structure is summarised in Fig.|| and, from (Q) and 
(0), in the small- wavenumber dispersion relation 



2lu = -(vi + V4:)ksin9± v a k ~ i(2D T + D L )k 2 
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between frequency to and wavenumber k, where 



v = y (vi — U4) 2 sin 2 9 + 4w 2 U3 cos 2 9, 
and 9 is the circular polar angle. 
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(19) 
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where A, B, C, and D are functions of a and a = |a| only. 
All dependence on the angle of tilt tf> is in the scalar and 
vector products in (|l7|). Evaluating the integrals, we find 
that B = = D (a consequence of the phase-only ap- 
proximation and the assumption a xy = 0), and A, C > 0. 



Inverting (17) we see that 



v u = M 



N ajaj 
ij ~ 1 + TV a 2 



(18) 
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where F = — — J f y is the Lorentz force, M — l/A, and 
N = C/A. (pf differs from that for a single vortex pP 



FIG. 2. The wave travelling along ±x that follows a local 
compression of an array of vortices moving along — y 
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We estimate the resulting wave-speeds for NbSe2 in the 
mixed phase jL2|. Its TDGL parameters are: \h ~ 700 
A, £ ~ 80l, T c w 7 K , = 5 ^ft cm. For an ap- 
plied transport current Jt = 1 A/cm 2 and inter- vortex 
separation a ~ A//, the wave-speeds c± ~ 1/zm/sec. 

The most obvious physical consequence of these waves 
is that the dynamic structure factor of a drifting flux lat- 
tice should display peaks at nonzero frequency (see (jl9|)). 
More dramatically, if a region of the flux-lattice moves 
past an impurity site, the impurity will "pluck" the flux 
lattice, and the effect will propagate along and trans- 
verse to the axis of drift, shaking up the lattice globally, 
through the sequence of events depicted in Fig. 2. This 
wave propagation in the absence of inertia is remarkable, 
and could well be a mechanism for nonthermal noise in 
drifting flux lattices. In addition, time-dependent exter- 
nal disturbances could excite resonances with the wave- 
like normal modes. 

Let us estimate the length scale l c above which 
these modes are actually propagative in character. For 
wavevectors k = (k x ,0) we see that 
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D L ~ MA [see after (Q), (g)] and v t ~ MF where 
F = Jt4>o/c is the Lorentz force per unit length on a 
vortex. Then 



A_ 
F' 



A » aH 2 /8ir [[13j, a being the flux-lattice spacing, so 
for a ~ 10/im and applied currents Jt ~ 100 A/cm 2 , 
£ c /a ~ 1, and the propagating modes should dominate. 
However, if a ~ 10 3 A, l c /a ~ 10 6 . 

In closing, we remark that our work settles an im- 
portant issue in the theory of the dynamics of moving 
flux lattices, namely their stability |14|. We have shown 
that dynamic interactions between vortices in a drifting 
flux lattice without inertia or pinning lead to a steady 
state with stable linear-response properties. Small dis- 
turbances about the drifting state travel as waves with 
a direction-dependent speed which, when calculated in 
terms of the parameters in the TDGL equations, turns 
out to be a few jum/s. These waves should be observ- 
able in systems with large flux-lattice spacing, at large 
imposed transport currents. The fast scanning tunneling 
microscopy approach of Troyanovskii et al. [15| seems to 
be the ideal way to observe these waves directly. 
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